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Abstract 



i In this paper, the semileptonic decays of heavy mesons are treated fully relativis- 

m ' 

On , 

| both at the inclusive and at the exclusive levels, and the spin of both parent and 

Oh, daughter particles are taken into account. The differential decay rates with respect 

Oh, to the lepton energy and momentum transfer are compared with data from ARGUS 

r-| . and CLEO. 



Introduction 

There are several reasons why semileptonic B decays are of interest. For one thing, the 
are a small variety of decay products, namely those which contain charm quarks (D, D*, 
D** , etc.) and those which do not (tt's etc.). V u b/V c b can be determined from the relative 
number of these decays. In addition, the heavy masses of the b and c quarks suggest 
that one might be able to apply perturbative QCD to calculate the strong corrections 
to these processes. This hope has been recently formalised in Heavy Quark Effective 
Theory Q. Finally, there's the fact that theoretical uncertainties are much smaller than in 
non-leptonic decays which contain a wider variety of hadronic decay products. 
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If we write out a parameterisation for the CKM Matrix, we see that it depends on a 
complex phase which is resposible for CP violation in the standard model. The magnitudes 
of the values for the CKM matrix elements place limits on the size of this phase and, thus, 
on the amount of CP violation in the standard model. According to the particle data 
group, V u b = .0035 ± .0015 and V c b = -040 ± .08. Recent values of V u b and V c b in the 
literature fall in this range [3-8] . 

In studying semi-leptonic decays, the first approximation is to neglect QCD and use a 
spectator model in which the up quark of the B meson is not involved in the decay except 
to recombine with the charm quark. In such a model, the decay of the B meson into a D 
meson reduces to that of the decay of a bottom quark into a charm quark. Quantities that 
can be determined directly from experimental data include the square of the momentum 
transfer, 

Q 2 = {B-Df 

= m 2 B + m 2 D -2E B E D + 2p B -p Dl (1) 

and the lepton energy, Ei. 

Inclusive Case 

Now that we have a process involving quanties that can be determined from experiment, 
we would want to come up with a theory that relates these quanties. One such model was 
the one devised by Altarelli, Cabbibo, Corbo, Maiani and Martinelli in 1982||. In this 
model, the bottom quark was assumed to be on shell and, thus, given by 

m 2 = (B-u) 2 

= m 2 B + m 2 u - 2m B \/pl+ml (2) 

in the B rest frame. This assumption has the advantage that it avoids having the decay rate 
depend on an arbitrary overall 1 / mjj that appears in a purely partonic treatment of these 
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decays fllPf The up quark momentum was then assumed to obey a Gaussian distribution 

1 (-P 
-^—exp — 

**pt V p f 



— 7— r f '•/' ( -Jr ) (3) 



which is normalised according to 

' d 3 P 0(p) = 1 (4) 



and is thus fixed up to an adjustable parameter, pf, known as the Fermi momentum. 

In our case, we want to consider the up quark as being more than a mere spectator: 
instead we write the effective uBb vertex as ^Vb{Pu)- The decay rate in this model is 

T{B - Udv) = 2m B (2^ J 2E C 2El 2E V 2K^ )|M| Sp-u-c-l-v) (5) 
where 

|Af| 2 = G 2 F L aP H al3 (6) 

and 

H ap = Tr [ 7 «(i _ l5 )(j, + mc ) 7 /3(l _ 7s )(£ + m b ) 75 (-^ + ro.hs^+rm)]. (7) 

It turns out this this integral is difficult to evaluate, the problem being that the up and 
charm quarks are not being assumed to combine into a specific meson. Instead, the quarks 
are combining to form a cluster X with four momentum X = u + c and mass m\ = (u+c) 2 . 
This mx is arbitrary save for the fact that, experimentally, mx > ttid — 1.8963 GeV 
while energy conservation requires that mx < m B - As a result, we will want to rewrite 
the hadronic phase space so that mx is integrated over this range. Using standard cluster 
decomposition techniques, the decay rate becomes: 

T(B -> XIV) = N]Vc ^ B? J \M\ 2 <t>(p u )-^d 4 Pc d 4 Pu d 4 (X u c)S(c 2 m 2 c )6(u 2 m 2 J 

L T d 4 Pl d 4 p u 6 4 (Q-l-v)6{l 2 )5(v 2 ) 



(2tt) 2 

-^d 4 Qd 4 XS 4 (B -Q- X)S(Q ■ Q - Q 2 )S(X 2 - m\ ) 
{2-kY 

dQ 2 dm 2 x 



(27T) 



3 



UsingjHj 



we get 



d 2 (X^ab) = (7r/2)X^(l,a 2 /X 2 ,b 2 /X 2 )^, (8) 



' d 4 Qd 4 XS 4 (B-Q-X)S(Q-Q-Q 2 )S(X 2 -m 2 x ) = (9) 



(2tt) 2 ^ v ^ / v-e / \ 27r2m B 

where pq = (1, X 2 /m 2 Bl Q 2 /m 2 B )mB = px- The remaining delta functions are 

5{c 2 -m 2 c ) = <5((V - U ) 2 - m 2 c ) 

= 6{X 2 +m 2 u -2E x E u + 2pxPuCos6 X u-m 2 c ) 

and 

= S((Q-lf) 

= S(Q 2 - 2E Q E l + 2p Q Eicos9 Q i) 

These cancel with the cosine integrations in d 3 p u and d 3 pi. The final expression for the 
decay rate is, thus, 



r ;> v;-, N \Vcb\ 2 V£ [ ^ PudpudEjdj) . 2 , ., 
T{B -> Xlv) = / \M\ <p{pu)—7. = rfQ rfm x (10) 



If we now compare this formula with data from ARGUS | TT| [ Y2\ and CLEO[^3| then we 
find, after minimising with respect to the parameters m u , m c , nib, Pf and |V^f,|/|14b|, that 
we get a good fit for parameters in the ranges 



m u 


= .13±.38GeV 


m c 


= 1.4 ± .4 GeV 


m b 


= 4.9 ± .3 GeV 


Pf 


= .5 ± .1 GeV 


and Vub/Vcb 


= .07 ±.05 



Note that the ARGUS and CLEO data include contributions from b — > civ and b — > 
ulV decays. In each case, the measured electrons were separated into different categories 
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including electrons from non-T(45) events, tp or ip(2S) decay, r decay or semileptonic D s 
decay, semileptonic D decay, 7r° — > e + e~ decay and semileptonic B decay, the latter being 
the ones that are used to make these plots. Additional background comes from having 
hadrons misidentified as electrons. Note that for Ei > 2.4 GeV, electrons from B decay 
can only come from charmless semi-leptonic decays. 

Using these parameters and taking the areas under the curves gives us the branching 
ratio Br(B -» cufi7)=10.09% and Br(B -» uulv) =.16%. Using 

T(B -> uclV) = Br(B -> uclv)/T B (11) 

and knowingQ that T B = (1.52 ± .11) x 10" 12 (1.52 x 10 24 ) GeV, \V cb \ can be calculated 
to be .034 ± .003. 

Exclusive Case 

In the spectator model, one can differentiate between D and D* mesons according to 
whether the daughter meson has spin or 1. That one can do this was overlooked in a 
recent paper by V. Barger et al that attempted to differentiate between different decay 
products in the differential mx distribution |p^. Mahiko Suzuki used this observation 
to calculate exclusive rates at zero Fermi momentum. In this frame, 

H ap = MqMq + MfMf (12) 

where 

A/ A oc Tr[^ + m)j x (l- l5 )(p + M)] /[4M{2m(E c + m)}i] (13) 

and 

Mi oc Tr[(t + m) 75 ^7 A (l- ia)(P + M)] /[4M{2m(E c + m)}?}. (14) 

Here e\ represents the three polarisations satisfying e\c x — 0. In the rest frame of c, is, 
therefore, given by 

e { p = (0,1,0,0), (0,0,1,0) 

e{ L) = (0,0,0,1). (15) 
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where (T) and (L) signify transverse and longitudinal polarisations, repectively. 

In the case where the b is not at rest in the B rest frame e is defined specifically in 
the B rest frame. Now, if instead of considering the light quark as a spectator, we treat 
it as an intermediate decay product in an effective theory involving a ubc loop, then the 
relevant traces are 

M A «Tr [^ + m c ) 7 A (l-7 5 )(|J + m ;) )75(-^ + m s )75] (16) 
M A ccTr [(t + m c yf x (l- l5 )(p + m b ) l5 (-i + m s )ft] (17) 

The Suzuki matrix elements are reproduced as p goes to zero. 
Starting with the B meson at rest, 

where 

N?G F V cb V B V D [ d 3 p b w/ , ,,1.. nn s 
S= ^ / 22gH0 (PuW(tu)\ 2 M (19) 

and where p u {f? u ) and are the up quark momenta in the B and D rest frames, respec- 

tively, Vb and Vd are the vertex constants and N is a normalisation. The wavefunctions 
<p(p u ) and tp(t u ) are 

4>{Pu) = -^exp[^A and = -^exp ( ^ j (20) 

where pf and tf are independent adjustable parameters. t u is given by 

P u = El -ml 

= [(E u E D -p u -p D )/m D ] 2 -ml (21) 

where E t is the energy of the up quark in the D rest frame. 
The phase space simplifies as follows: 

d 3 (B - DW) cx ^^M^^^SHB -D-l-u) (22) 

- f fefeSj dcose^s^dp^ldcose^dp^dcose'^ (23) 
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B — D — I — v = and v 2 



implies 



(B - D - I) 2 = 
=> 2p D Eicos<t>t = 2E D E l - 2m B (E D + E t ) + {m 2 D + m%) 



(24) 



(25) 



so the phase space becomes 



d s (B -» Dlu) oc 



it 2 Pd dp D dEi 



dp u pl dcos6 u dtp u dp' u p' 2 dcos9' u dtp' u 



(26) 



8Ed E u E' u 



If we now insert into this model the parameters given in the previous section we run into 
problems: it turns out that we only get agreement with ARGUS |Q data for low values 
of Q 2 . This is presumably due to final-state interactions, which in a perturbative QCD 
framework are expected to grow as one approaches the end-point of the Q 2 distribution. In 
the exclusive case, QCD corrections are restricted to those which do not create additional 
hadrons, that is quark propagator self-corrections and vertex corrections. Corrections 
to the cDu vertex are of particular interest because they provide a phenomenological 
explanation for the discrepancy: the exchange of a gluon between the up and charm quark 
can reduce their relative momentum, allowing them to combine to form a D or D* meson. 

Conclusion 

This model effectively describes the dependence of both inclusive and exclusive semilep- 
tonic B decays on the Fermi momentum of the constituent quarks. The parameters that 
arise naturally in this model agree with those used in other models. 

Given that this model describes both inclusive and exclusive decays, we can estimate 
the rate of semileptonic B decays into D** mesons or clusters consisting of Z?'s or D* 
and 7r's by subtracting the exclusive rates into D and D* from the inclusive semileptonic 
B decays into charmed mesons. Experimentally, this rate is found to be between 33% 
and 41% of the total semileptonic rate|l^|(l^]. This model would appear to have the best 
chance of accounting for all possible semileptonic decay products of B mesons. 



Figure 1: |§P for m»=.13 GeV, m c =l.A GeV, m 6 =4.9 GeV, p f =.5 GeV and V ub /V cb «.07 
with data from ARGUS [11] [12] and CLEO[13] 
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Figure 2: for m„=.13 GeV, m c =1.4 GeV, ™ fc =4.9 GeV,p/=.5 GeV and V ub /V cb «.07 
with data from ARGUS [17] 

References 

[1] Particle Data Group, Phys. Rev. D50 (1994), 1173. 

[2] M. Neubcrt, CERN Preprint #7396 (1994). 

[3] H.-N. Li, and H.-L. Yu, Chung-Cheng University Preprint (1994). 

[4] M. Neubert, CERN Preprint #7396 (1994). 

[5] V. V. Kisclcv, IHEP Preprint #94-77 (1994). 

[6] M. G. Olsson and S. Veseli, Madison Physics Preprint #851. 

[7] A. A. El-Hady, K. S. Gupta, A. J. Sommerer, J. Spence and J. P. Vary, Iowa State 
University Preprint (1994). 

[8] C. S. Kim and A. D. Martin, Proceedings of the International Workshop on B-Physics, 
Nagoya, Japan (1994). 

[9] C. Altarelli, N. Cabbibo, G. Corbo, L. Maiani and G. Martinclli, Nuclear Physics 
B208 (1982), 365. 

[10] V. D. Barger and R. J. N. Phillips, Collider Physics, Addison- Wesley Publishing 
Company, California (1987). 

[11] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B249 (1990), 359-365. 

[12] ARGUS Collaboration, H. Albrecht ct al, Phys. Lett. B318 (1993), 397-404. 

[13] CLEO Collaboration, P. Avery et al CLEO Conference #94-7 (1994). 

[14] V. Barger, C. S. Kim and R. J. N. Phillips, Madison Physics Preprint #501 (1989) 

[15] M. Suzuki, Nuclear Physics B258 (1985), 553. 

8 



[16] ARGUS Collaboration, H. Albrecht et al, Z. Phys. C57 (1993), 533-540. 
[17] ALEPH Collaboration, D. Buskulic et al, CERN Preprint #94-173 (1994). 
[18] OPAL Collaboration, R. Akcrs et al, CERN Preprint #95-02 (1995). 



9 



